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Blasius boundary layer evolution is studied by means of bicoherence calculations. The layer 
is acoustically excited at the T-S frequency to provide a controlled transition. Measurements 
are made using a smooth surface as well as various roughness patterns. The bicoherence 
calculations are used to determine the extent to which frequency resonant velocity 
fluctuation waves can participate in energy exchange. The emphasis is on downstream 
variation of the individual interactions among harmonic modes. A limited picture of the role 
of quadratic wave interactions is revealed. 
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digital Fourier transform of the streamwise velocity fluctuation 
roughness element thickness 
acoustic forcing level 

roughness element wavelength in the x direction 

kinematic viscosity 

roughness element orientation angle 


I. Introduction 

S kin friction and heat transfer rates are significantly higher in turbulent flow. The transition to turbulence occurs 
upon the onset of non-linear interactions. Advanced time series techniques such as the calculation of the 
bicoherence are conceived to determine the nature of quadratic interactions. In the early transitioning flow the 
fundamental Tollmien-Schlichting wave distorts and generates harmonics that interact with one another. King and 
Breuer 1 discuss the mechanism of receptivity by which external disturbances determine the wavelike disturbances in 
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the early flow and show amplitude growth data for various acoustic disturbance levels and wall roughness patterns. 
Results obtained by Davila, King and Chokani 2 reveal that the early non-linear interactions persist downstream in 
narrow frequency ranges. This discreteness makes it difficult to present the data using the often-shown contour 
graphs. This is especially true for smooth surfaces. In the present work several smooth- and rough-walled cases are 
studied for acoustically excited transitioning layers by examining the bicoherence value as a function of downstream 
position for selected interactions. This represents a departure from the more common presentation, away from the 
panoramic bispectral picture, in favor of individual wave triads. 

II. Experimental Set-Up 

The experiments were conducted in the closed-loop, low-turbulence (2 ft. by 3 ft.) Low-Speed Wind Tunnel 
located at NASA Langley Research Center at a freestream velocity of U x ~ 11.1 m/s. The wind tunnel has a 6.1 m 
long test section with a 0.61 m (high) by 0.91 m (wide) cross-section. The contraction ratio is 10:1 and the 
measured turbulence intensities are approximately 0.1% in the range of 0. 1 < f < 400 Hz. The boundary layer under 
investigation developed on a 2.62 m long, 12.7 mm thick, aluminum plate, with a 24:1 elliptical leading edge, which 
spanned the width of the test section. The measurement surface was polished to 0.2 pm r.m.s. surface finish. A 
trailing edge flap was used to position the stagnation point on the measurement surface. 

A rectangular cutout, located 24.6 cm from the leading edge of the plate, is equipped with interchangeable 
inserts to allow mounting surface “roughness” pattern about Branch I. The patterns consist of a series of thin 
rectangular strips attached to the insert surface. In all cases, they were of equal width, spaced by their own width, 
and having a constant thickness Ah (elevation above the plate surface); this produces a square -wave roughness 
pattern. They were generated using a photolithographic process. Only the medium and high roughness cases, with 
Ah = 36 pm, and Ah = 72 pm, respectively, are shown. Results were also obtained with Ah = 18 pm, but they 
showed no special features. The angle ip is the orientation of the length of the strips with respect to the y axis (the 
spanwise direction). It is equivalently the misalignment of the surface roughness wavenumber vector with respect to 
the x axis (the streamwise direction). In patterns with vp — 0° the surface roughness wavenumber points in the 
direction of flow, thus being two-dimensional. In patterns with \p > 0° the wavenumber vector has a component in 
the spanwise direction; thus a three-dimensionality is introduced. The x-component of the wavenumber has an 
associated x-direction wavelength /, x that defines the spatial periodicity of the pattern in the streamwise direction. 
Table 1 summarizes the geometric details of the patterns. The values of /, x were chosen to be close to the T-S 
wavelength, 50.2 mm at Branch I. Details of the wind tunnel and the roughness patterns are given in reference 1. 


pattern 

X x [mm] 

V 

Ah [pm] 

C (2-dimensional) 

50.3 

0° 

72 

B (2-dimensional) 

50.3 

0° 

36 

M (3-dimensional) 

50.7 

15° 

36 

K (3-dimensional) 

52.4 

30° 

36 

N (3-dimensional) 

55.5 

45° 

36 


Table 1. Surface roughness pattern geometry details. 


Five 203.3 mm diameter woofers were used for acoustic excitation; four speakers were mounted on the upstream 
wall and the other on the downstream wall. The amplitude of the acoustic field, u ac , could be varied to produce 
different levels of acoustic forcing, s = u /U^. More complete details of the experiment are given by King, 3 and a 
select group of acoustic and roughness combinations are presented here. Low acoustic levels produced no 
observable bicoherence. Only medium and high excitation levels, with s = 3.8 x 10' 5 and s = 7.6 x 10" 5 , 
respectively, are discussed here. 

The phase-locked velocity fluctuations were derived from the signal of a single-element hot-wire, operated at an 
overheat ratio of 1.8 with a constant temperature anemometer. The acquisition of the hot-wire signal was triggered 
with the TTL of a reference signal and sampled at a frequency of 1 kHz. The reference signal was the excitation, or 
forcing frequency, f 0 , which was applied to the loudspeakers that were used to generate the downstream traveling 
acoustic wave. The dimensionless frequency of the T-S wave was F = 55 x 10' (1 (=2jrfv/Uaf), which translates to a 
physical forcing frequency of f„ = 71 Hz. Results are shown for four streamwise stations, given by modified 
Reynolds number R (= \ : Rc x ) = 943, 1038, 1124 and 1205 with the smooth and rough plates, two-dimensional and 
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three-dimensional. While the data are limited, an attempt is made to discern features of the quadratic nonlinearities 
that are determined by acoustic excitation level and by roughness element thickness and orientation. 

III. Calculations 

The nonlinear energy cascade that occurs due to the three-wave resonant interactions was quantified using the 
bicoherence. Kim and Powers 4 give a detailed description of the calculation technique for the bicoherence. Corke 
and Mangano 5 , Miksad et a 1 . ' among others have used the bicoherence to study transition to turbulence, but a few 
salient details are repeated here for the sake of completeness. The bispectrum is a third-order cumulative function 
that is defined as: 

S b (f i ’ f j ) = e[ X * (f ; )X * (f . )X(f ; + f j )] (1) 


where E[ ■] is the expected value, X(f) are the digital Fourier transforms of the time series records, the superscript 
asterisk denotes the complex conjugate. The normalized representation of the bispectrum is the bicoherence, which 
is bound by 0 and 1. The bicoherence square is defined by: 
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If the triad that is present at f. , f. and the sum mode f.+f. is comprised of independently generated oscillations, 
then these are statistically independent of each other, and consequently b 2 (f, f) ~ 0.0. On the other hand, high 
values of the bicoherence indicate a high degree of temporal phase coherence between the product of the oscillations 
(at f. and f ) and the oscillation at the sum frequency (at f+f). The sum interactions are of the form 

(f. , f) -► f+f. 

v i’ y i j 

The symbol — * denotes “phase -coupled interaction with the mode at.” The targeted mode is determined by adding 
the ordinate frequency to the abscissa frequency in a triangular region defined by the vertices (0, 0), (f n , 0), 
(f/2, f/2), where f is the Nyquist frequency, 500 Hz. This is the top triangular zone of the plot. The difference 
interactions are of the form, 

(f. , -f.) -> f.-f 

which are the ordinate frequency subtracted from the abscissa frequency, are in the region defined by the vertices 
(0, 0), (f n , 0), (f n , -f n ). This is the bottom triangular zone of the plot. The bicoherence square, b 2 (f , f), with sum and 
difference interactions is often plotted as a contour level in the (f., f. ) plane. Figure 1 (reproduced from reference 2) 
is representative of such plots; it is shown here for illustration. 

Figure 1 is not typical. In the present results, most contour plots showed narrow peaks involving mostly the 
fundamental and the harmonics. These are not easily observable; therefore, on this occasion only ten triads, all 
involving harmonic peaks, were selected for observation. The interactions are defined using the relative frequency 
f/f 0 , where f is the frequency and f 0 is the fundamental frequency of f Q = 71 Hz. The fundamental 71 Hz and the 
harmonics 142 Hz, 213 Hz and 284 Hz correspond to f/f„ values of 1, 2, 3 and 4. For convenience, in the discussion 
of results and in all figures, the bicoherence square is denoted by B, in place of b 2 . For example, B(l,2) = b 2 (f 0 , 2f 0 ), 
corresponding to a sum interaction targeting the frequency 3f 0 . Similarly, B(3, -1) = b 2 (3f G , -f Q ), corresponding to a 
difference interaction targeting the frequency 2f 0 . 

The number of downstream positions where measurements were made precludes proper observation of the 
expected widening of the bicoherence contours as turbulence breakdown sets in. But in fact, turbulence breakdown 
does not occur in most of these cases. Instead, the present discussion will focus the downstream variation of the 
bicoherence values for the selected triads and the corresponding evolution of the energy spectrum. 

The boundary layer is excited using an unstable mode of 71 Hz. This forced fundamental mode is clearly 
identifiable at all streamwise locations; the 60-Hz signal (and its harmonics), due to signals picked up from the 
electric cables, can also be seen as a series of narrow spikes. 
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Figure 1. Bicoherence square (top) and Energy Spectrum at R=1082 for the 
case of high acoustic exitation and high roughness, pattern C. 
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Case: Medium excitation, smooth plate 
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Figure 2. Bicoherence and energy spectrum as a function of R for the medium-excitation, smooth-plate case. 


IV. Results 

Results for for the case of medium acoustic excitation and smooth plate case are shown in figure 2. The energy 
spectrum is uneventful, and the fundamental and harmonic mode amplitudes are low. The low amplitude of the 
fundamental prevents it from transferring sufficient energy to the harmonic modes. All triads show increasing 
bicoherence values up to R = 1124, but the modes at f/f Q = 2, 3, and 4 show low energy values at all stations, 
indicating poor energy transfer. As a result, at R = 1205 the harmonic modes randomize, resulting in lower 
bicoherence values for all triads. With low bicoherence values at R = 1205, energy transfer is inhibited, and 
therefore the transition does not occur in the R region studied. 
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Case: High excitation, smooth plate 
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Figure 3. Bicoherence and energy spectrum as a function of R for the high-excitation, smooth-plate case. 


Figure 3 shows the results for the case of high acoustic excitation with a smooth plate. As in the previous case, 
an uneventful energy spectrum shows that, without surface roughness, the amplitude of the fundamental mode 
remains low at all stations, despite the high excitation. The triads (1,1) and (2,-1) show high bicoherence due to the 
growth of the mode at f/f 0 = 2. All other triads reach bicoherence values near 0.9. However, energy transfer to the 
harmonic modes is not effective, as evidenced by their persistently low harmonic energy levels. With low amplitude 
levels in the harmonic modes, the bicoherence has decreased for all triads at R = 1205. The transition does not occur 
in the R region studied. In fact, in both smooth-plate cases, the acoustic excitation alone is not enough to induce 
transition in the R region studied. 
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Case: Medium excitation, high roughness, \j/ = 0° (pattern C) 
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Figure 4. Bicoherence and energy spectrum as a function of R for the medium-excitation, high-roughness, 
\|/ = 0° (pattern C). 

The results for the medium excitation, high roughness case are shown in figure 4. Evidently the introduction of 
surface roughness in the receptivity region has resulted in an increase of the fundamental mode amplitude 
downstream. This behavior is consistent with measurements by King and Breuer 1 for two-dimensional roughness. 
The harmonic mode amplitudes are also greater. At R = 1124 all triads have high bicoherence, and the harmonic 
modes (up to f7f 0 = 4) have grown more than they did in either smooth-plate case, which is evidence of effective 
energy transfer. At R = 1205, the filling of the energy valleys shows that transition is incipient. As a result, the 
higher-frequency harmonic modes begin to randomize, and the bicoherence values begin to decrease. At the last 
station, for all triads that include the mode at f/f 0 = 4, the bicoherence value falls below 0.4. 
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Case: High excitation, high roughness, \j/ = 0° (pattern C) 
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Figure 5. Bicoherence and energy spectrum as a function of R for the high-excitation, high-roughness, \j/ = 0° 
case (pattern C). 


The results for the high-excitation, high-roughness case are shown in figure 5. At all stations the fundamental 
and harmonic modes have higher amplitude than in the high-excitation, smooth case. In addition, all triads reach 
high values at R = 1038, indicating early strong interactions with effective energy transfer to harmonic modes up to 
f/f 0 = 4. Transition has begun at R = 1 124 and the layer is turbulent at R = 1205. The effect of the roughness region 
is evident: increased mode growth, stronger interactions, effective energy transfer, and early transition. 

Data for many patterns were collected at medium acoustic excitation levels. These are shown in the remaining 
results. 
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Case: Medium excitation, medium roughness, \j/ = 0° (pattern B) 
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Figure 6. Bicoherence and energy spectrum as a function of R for the medium-excitation, medium-roughness 
case, \|/ = 0° (pattern B). 

The results for the medium-excitation, medium-roughness case are shown in figure 6. Relative to the high- 
excitation, smooth-plate case (Fig. 3), the fundamental-mode amplitude is greater. The amplitudes of the modes at 
f/f 0 = 2 and 3 are also greater, suggesting more effective energy transfer. In fact, for the (1,1) and (2,-1) triads, high 
bicoherence values persist farther, despite the lower excitation level. Roughness is evidently determinant here. 

Relative to the medium-excitation, high-roughness case (Fig. 4), the medium-excitation, medium-roughness case 
shows lower fundamental and harmonic mode amplitudes, as well as lower bicoherence values for all but triads (1,1) 
and (2,-1). The mode at f/f 0 = 3 loses amplitude, and the mode at f/f 0 = 4 does not receive energy. Evidently there is 
no effective energy transfer to the higher harmonics, and therefore transition does not occur in the R region studied. 
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Case: Medium excitation, medium roughness, i j/ = 15° (pattern M) 




Figure 7. Bicoherence and energy spectrum as a function of R for the medium-excitation, medium-roughness 
case v|/ = 15° (pattern M). 

Figures 7, 8 and 9 show the results of the three-dimensional cases (obtained using with oblique surface roughness 
patterns) with \|/ = 15°, 30 0 and 45°, respectively. The energy spectrum plots show that the fundamental mode 
amplitude decreases with increasing i|/. This behavior is consistent with the measurements by King and Breuer 1 for 
oblique surface roughness. The case with \|/ = 15° (Fig. 7) is similar to the two-dimensional case of figure 6. All but 
the (1,1) and (2,-1) triads have low bicoherence values. For the case with \|/ = 30° (Fig. 8) the bicoherence values of 
triads (1,1) and (2,-1) are lower at R = 1 124 and 1205. The effect is yet more pronounced for the case with \|/ = 45°, 
(Fig. 9) with low bicoherence of triads (1,1) and (2,-1) at all but the first station. 
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Case: Medium excitation, medium roughness, \| / = 30° (pattern K) 




Figure 8. Bicoherence and energy spectrum as a function of R for the medium-excitation, medium-roughness 
case, \|/ = 30° (pattern K). 

In general, the bicoherence results reveal details of modal evolution not evident just from the energy spectrum. 
Two-dimensional roughness yields increased fundamental mode amplitude, and facilitates quadratic interactions that 
can lead to transition. To wit, the one case (Fig. 4) in which transition to turbulence appears incipient and the one 
case (Fig. 5) in which it occurs, are of high roughness. In both, the transition is preceded by strong interactions of 
all triads. In the medium excitation, medium roughness cases (Figs. 6 to 9) the energy spectrum and bicoherence 
results show that a larger orientation angle \|/ inhibits fundamental-mode growth and quadratic interaction of the 
selected triads. 
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Case: Medium excitation, medium roughness, \|/ = 45° (pattern N) 




Figure 9. Bicoherence and energy spectrum as a function of R for the medium-excitation, medium-roughness 
case, v| / = 45° (pattern N). 


V. Conclusion 

The results show that the roughness element thickness and orientation influence the development of quadratic 
interactions. Only part of the picture is evident from the evolution of energy spectrum values. In conjunction with 
these, bicoherence calculations can reveal a more complete picture of the evolution of interacting triads. 
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In particular, the energy spectrum and bicoherence calculations show that: (1) in the two-dimensional case, 
higher surface roughness yields higher fundamental mode amplitudes and promotes quadratic interactions, and (2) 
that three-dimensionality imposed by surface roughness obliqueness inhibits quadratic interactions. The results 
confirm those of King and Breuer 1 . Also, they suggest that the presence of multiple, effective quadratic interactions 
are a precondition for transition in the R region studied. 
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